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ABSTRACT  (GmHmm  m  CM  If  wmimr  M  I9mtttr  ftp  ftfMfc  wiif) 

A  new  satellite  triaxial  accelerometer  system  has  been  developed.  This 
Instrument  has  been  flight-tested  on  three -axis  stabilized  satellites  and  has 
demonstrated  the  capability  to  accurately  measure  accelerations  in  the  satel¬ 
lite’s  in-track,  cross -track,  and  radial  directions.  The  in -track  data  provide 
direct  determination  of  atmospheric  density.  These  data,  obtained  during  a 
period  of  high  solar  flux,  supplement  the  extensive  set  of  measurements  ob¬ 
tained  by  other  accelerometer  experiments  during  low  solar  flux  conditions. 
The  cross-track  and  radial  results  permit,  for  the  first  time,  large-scale 


20.  Abstract  (Continued) 

^measurements  of  the  zonal  and,  possibly,  radial  components  of  neutral 
atmospheric  winds.  These  simultaneous  wind  and  density  data  provide  an 
extremely  valuable  input  for  the  understanding  of  dynamic  processes  in  die 
atmosphere  and  for  the  improvement  of  atmospheric  models. 
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Preface 


Funds  to  Initiate  the  accelerometer  development  were  provided  by  the  AFGL 
In-House  Independent  Laboratory  Director's  Fund.  Murray  Meldrum  and  William 
Lange,  Bell  Aerospace  Textron,  were  principal  investigators  for  instrument  fabri¬ 
cation  under  AF  Contract  F10828-74-C-O114. 

At  NS WC  Linda  Lynch  and  Becky  Boyer  were  responsible  for  daily  quick-look 
data  analysis  and  generation  of  raw  data  tapes  for  AFGL.  Additional  assistance  in 
data  quality  analysis  was  provided  by  Carol  Brandi.  The  data  processing  and 
plotting  software  was  generated  by  Allen  Fisher,  Ling  Szeto,  and  Dee  Batayte. 

Processing  of  density  and  wind  data  was  accomplished  by  R.  W.  Fioretti,  A.  J. 
Mazzella,  Jr. ,  and  R.  Leong  of  RDP,  Inc.  under  AF  Contract  F19628-76-C-0244. 
Material  for  Figures  14  and  17b  was  generated  by  Edward  Robinson  (SUWA). 

Dorothy  Gillette  and  Sammy  Chan  assisted  in  various  aspects  of  data  presentation. 
Margaret  Anderson  and  Stephanie  Fur  low  typed  the  numerous  drafts  and  final  manu¬ 
script.  Critical  technical  comments  were  provided  by  K.S.W.  Champion  and  J.  M. 
Forbes. 
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Application  of  the  Satellite  Triaxial 
Accelerometer  Experiment  to  Atmospheric 
Density  and  Wind  Studies 


1.  INTRODUCTION 


A  single-axis  accelerometer  system,  MESA  (Miniature  Electrostatic  Acceler¬ 
ometer),  has  been  flown  successfully  by  the  Air  Force  Geophysics  Laboratory 
(AFGL)  on  eight  satellites.  Table  1  lists  the  spacecraft,  their  launch  dates  and 
lifetimes.  Except  for  S3-1  and  S3-4,  all  flights  utilized  three  single-axis  sensors 
mounted  orthogonally  to  achieve  a  triaxial  sensing  capability.  The  resulting  pack- 
age  typically  weighed  10  kg  and  consumed  about  20  W  of  power.  The  accelerometer 
proved  to  be  extremely  accurate  and  reliable  on  all  missions.  For  example,  on 
the  Atmosphere  Explorer -C  (AE-C)  satellite,  accelerations  were  measured  over 
a  five-year  period.  During  its  first  year  of  operation,  with  the  spacecraft  in  a 
highly  elliptical  orbit,  the  Instrument  was  used  to  measure  density,  orbit-adjust 
thrusts,  and  vehicle  motions.  For  about  the  next  four  years  AE-C  was  in  a  circu¬ 
lar  orbit  too  high  for  accurate  density  data  to  be  obtained.  However,  the  instru¬ 
ment  was  operated  as  required  for  monitoring  of  thrusts  and  vehicle  motions  by 
NASA  to  control  the  satellite.  During  the  last  few  weeks  of  the  AE-C  lifetime, 
when  its  altitude  was  below  250  km,  density  measurements  were  again  obtained. 

The  AE-E  triaxial  accelerometer,  launched  in  1975,  has  been  utilized  in  a  similar 
manner,  and  reentry  data  were  obtained  during  June  1981.  Marcos  et  al1  have 


(Received  for  publication  2  March  1982) 

1.  Marcos,  F.A. ,  Mclnerney,  R.E.,  and  Fioretti,  R.W.  (1978)  Variability  of  the 
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A  complete  desdnrt&n  of  MESA  Instrument  that  was 


flown  successfully  on  three  Atmosphere  Explor<  r  satellites  has  been  given  by 
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Champion  and  Marcos.  To  develop  a  triaxial  sensing  capability  the  suspension 
system  was  modified  by  instrumenting  both  radial  axes  with  precision  constrain- 
ment  loops.  A  brief  description  is  given  in  the  following  paragraphs. 

Figure  1  shows  die  mechanical  assembly  of  die  instrument.  The  one  gram 
proof  mass,  in  the  center  of  the  assembly,  is  constrained  in  all  three  axes  by 
means  of  electrostatic  forces  generated  by  electrodes.  There  are  a  total  of  14 
electrodes,  three  on  each  of  the  forcer  rings  facing  both  sides  of  the  flange  and 
eight  on  the  carrier  cylinder  facing  the  inside  of  the  proof-mass  cylinder.  The 
Z  axis  is  along  the  axial  centerline  of  the  proof-mass  cylinder  as  in  the  single-axis 
instrument.  The  X  and  Y  axes  are  perpendicular  to  the  cylinder  and  are  each 
defined  by  the  four  pairs  of  electrodes  on  the  carrier.  The  proof -mass  cylinder 


;  .  .... 


Figure  1.  Instrument  Mechanical  Assembly 


2.  Champion,  K.S.W.,  and  Marcos,  F,A.  (1973)  The  triaxial  accelerometer 
system  on  Atmosphere  Explorer,  Radio  Sci.  8: 197. 
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is  0.  025-cm  thick  and  its  flange  is  0.  0375-cm  thick.  The  spacer  has  a  thickness 
of  0.  056  cm.  The  gap  between  the  flange  and  each  Z-axis  forcer  electrode  is 
0.  0089  cm.  Between  the  X-  and  Y-axes  carrier  electrodes  and  the  inside  of  the 
proof -mass  cylinder  the  gap  is  0.  0033  cm. 

Accelerations  along  the  Z  axis  are  measured  by  the  set  of  six  electrodes  fac¬ 
ing  the  flange.  The  X  and  Y  axes  are  each  controlled  by  two  pairs  of  electrodes  on 
opposite  sides  of  the  carrier.  The  X  and  Y  axes  constrainment  systems  are  identi¬ 
cal,  each  using  two  constrainment  circuits  to  prevent  rotation  in  two  axes.  The  Z 
axis  requires  only  a  single  circuit.  The  proof  mass  is  not  constrained  against 
rotation  about  the  cylinder  axis.  Figures  2  and  3  are  block  diagrams,  for  the  Y 
(or  X)  and  Z  axes  respectively,  of  the  constrainment  loop  connections  to  the  elec¬ 
trodes  within  the  mechanical  assembly.  Measurement  of  acceleration  is  accom¬ 
plished  by  sensing  motion  of  the  proof  mass  with  respect  to  the  instrument  case. 

An  electrostatic  force  is  then  generated  on  the  appropriate  electrodes  to  restore  the 
proof  mass  to  its  null  or  zero  input  position.  The  restoring  force  is  generated  by 
a  dc  voltage  proportional  to  input  acceleration.  This  voltage  is  fed  to  a  precision 
A/D  converter.  The  output  of  the  instrument  consists  of  a  digital  pulse  rate  pro¬ 
portional  to  the  input  acceleration. 


Figure  2.  Block  Diagram  of  Y-  (or  X-)  Axis  Force  Rebalance  Loop 


£-ii 


The  magnitude  of  the  accelerations  that  must  be  measured  varies  from 
-2 

1.2  X  10  g  for  a  maximum  orbit-adjust  engine  thrust  down  to  approximately 

.a 

1  X  10  g  for  high -altitude  drag  data.  Three  sensitivity  ranges  are  used  for  each 
axis  to  cover  this  measurement  interval.  These  ranges  are  selected  automatically 
during  flight.  The  scaling  is  selected  independently  for  each  axis.  The  nominal 
maximum  g  level  for  each  const rainment  range  is  given  in  Table  2  for  a  test  flight 
unit  and  the  operational  flight  unit. 

The  instrument  maximum  pulse  rate  output  representing  full  scale  in  a  posi¬ 
tive  or  negative  direction  is  IS  kHz.  The  digital  accumulator  is  a  bidirectional 
or  up-down  counter  that  stores  the  incoming  pulses  for  a  given  sampling  time  and 
then  provides  three  output  words  in  serial  form.  ThiB  output  word  represents  die 
average  acceleration  input  over  the  2  .  045-sec  sample  period.  Instrument  temper¬ 
ature  is  monitored  as  an  analog  word  once  every  114.  5  sec  by  a  thermistor  placed 
near  the  proof  mass.  Based  upon  this  pulse  rate  and  sample  time,  accelerations 
as  low  as  2X  10”9g  can  be  sensed  on  the  most  sensitive  scale  of  the  operational 
unit. 
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Table  2.  Accelerometer  Sensitivity  Ranges 


Test/ Future 
Flights 

Operational 

Flight 

15000  pg 

6000  pg 

1500 

600 

300 

60 

15000  Jig 

15000  Mg 

1500 

900 

150 

60 

Z  Axis 
Ranges  A 
B 


An  analysis  of  satellite  mechanical  vibrations  Induced  by  the  attitude  control 
3  4 

thrusters  was  performed.  *  These  studies  indicated  noise  peaks  at  several  fre¬ 
quencies.  Electronic  filtering  was  added  to  the  accelerometer  to  minimize  die 
effects  of  the  noise  accelerations.  The  bandwidth  response  is  as  follows:  1. 0  at 
dc,  0.  997  at  0. 1  Hz,  0.  5  at  0. 2 5  Hz,  and  0.  003  at  6.  0  Hz. 

Accelerometer  scale  factor  calibration  is  achieved  by  introducing  a  fraction 
of  the  earth's  lg  field  into  the  sensitive  axis.  This  is  accomplished  by  tilting  the 
sensitive  axis  by  a  known  angle  from  the  horizontal,  zero  input  position.  Accuracy 
of  the  scale  factor  measurement  is  limited  by  the  ability  to  generate  and  maintain 
very  small  angular  displacements.  To  achieve  ground  calibration  data  at  very  low 

input  levels  a  laser  interferometer  is  used.  It  has  a  resolution  capability  of 
•7 

0.  02  arc  sec  (10  g).  The  scale  factor  is  a  function  of  temperature  due  mainly 
to  the  expansion  and  contraction  of  mechanical  parts.  Temperature  calibrations 
were  performed  in  the  range  between  70°  and  120°F.  The  calibration  involved 
repeated  measurements  with  emphasis  on  the  Z  axis  and  its  most  sensitive  scale 
(C-range  of  Table  2).  Standard  deviations  of  the  C  range  calibration  data  were 
typically  0. 1  Mg.  Scale  factor  temperature  coefficients  were  of  the  order  of 
0. 1  percent/°F. 

The  accelerometer  flight  hardware  is  shown  in  Figure  4.  The  instrument 
weltfis  4.  6  kg,  and  has  the  dimensions  8.  9  X  12.7  X  29.2  cm.  Maximum  power 

3.  Counts,  J.  (1975)  Effects  of  Flexible  Solar  Arrays  on  Sensed  Accelerations  in 

a  Rigid  Spacecraft.  VPI-E-75-13,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  Virginia. 

4.  Butler,  T.A.  (1977)  Effect  of  Attitude  Control  Thruster -Induced  Structural 


Vibrations  on  Sensed  Accelerations  of  a  Space  Vehicle.  Los  Alamos 
boratories,  Informal  Report  LA -6965-MS,  Los  Alamos, 
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Figure  4.  Accelerometer  Flight  Package 


consumption  is  11. 7  W.  The  accelerometer  is  mounted  with  its  Z  axis  aligned 
with  the  satellite  in-track  direction,  the  Y  axis  along  the  earth's  radial  direction, 
and  the  X  axis  in  the  cross-track  direction.  Due  to  satellite  physical  limitations 
the  instrument  was  mounted  below,  to  one  side,  and  forward  of  the  center  of  gravity. 

The  telemetry  data  tapes  were  processed  at  NSWC  for  de commutation  and  time- 
tagging.  Tapes  containing  only  time,  accelerometer  raw  data,  and  temperature 
were  then  prepared  and  shipped  routinely  to  AFGL  for  processing  to  derive  density 
and  wind  data.  Experiment  operation  was  monitored  at  NSWC  by  preparing  plots 
of  acceleration  vs  time  on  a  routine  basis. 


3.  ORBITAL  FLIGHT  DATA 

An  orbital  accelerometer  flight  test  was  conducted.  After  turn-on  die  auto¬ 
matic  range  circuit  functioned  normally  and  put  all  three  axes  into  their  most 
sensitive  range.  This  indicated  that  the  preflight  estimates  of  the  orbital  accel¬ 
eration  environment  and  instrument  response  characteristics  were  accurate.  Data 
obtained  on  all  three  axes  during  one  orbit  are  shown  in  Figure  5.  The  top  frame 
stows  the  accelerations  sensed  in  the  along-track  (Z)  direction.  The  next  two 
frames  show  accelerations  measured  in  the  radial  (Y)  and  cross-track  (X)  direc¬ 
tions,  respectively.  The  bottom  frame  shows  the  temperature  (°F)  measured 
inside  die  proof-mass  housing.  A  typical  curve  of  the  Z-axis  drag  acceleration  vs 
time  would  vary  at  the  orbital  rate  and  have  its  maximum  value  at  the  time  die 


satellite  is  at  its  lowest  altitude  (where  the  density  is  highest).  Drag  accelerations 
are  sensed  with  a  negative  polarity  for  tills  experiment.  Hence,  maximum  densi¬ 
ties  correspond  to  more  negative  acceleration  values  and  minimum  densities  cor¬ 
respond  to  less  negative  acceleration  values.  For  the  data  in  Figure  S.  the  instru¬ 
ment  bias  shifts  the  curve  downward  for  the  Z  axis.  The  instrument  successfully 
provided  measurements  on  all  three  axes.  Several  anomalies  were  observed  and 
are  discussed  below. 

Examination  of  the  in-track  data  in  Figure  5  shows  general  characteristics  of 
a  drag  profile  with  the  largest  drag  near  minimum  altitude.  However  the  results 
are  dominated  by  an  opposite  polarity  spike  peaking  at  65400  sec  GMT.  Hits 
spike  has  a  2-pg  amplitude  and  about  500-sec  duration.  It  was  subsequently  de¬ 
termined  that  this  most  significant  deviation  from  the  anticipated  drag  profile  was 
caused  by  the  large  heat  input  experienced  when  the  instrument  was  exposed  to 
solar  radiation.  Onset  occurred  when  the  satellite  exited  penumbra.  At  this  time 
the  instrument  temperature  begins  to  increase.  The  result  of  tills  temperature 
increase  is  to  cause  an  apparent  decrease  in  sensed  acceleration.  Termination 
of  the  spike  occurred  when  the  vehicle  pitched  down  sufficiently  to  shade  the  in¬ 
strument.  Further  study  revealed  that  the  accelerometer  was  responding  in  a 
complicated  manner  to  thermal  variations  throughout  the  orbit.  Other  anomalies 
found  in  the  data  were  the  following: 

(a)  Bias  values  were  initially  high  on  all  three  axes.  The  data  in  Figure  5 
shows  a  bias  of  approximately  -20  fig,  -2  fig,  and  +12  fig  for  the  Z,  Y,  and  X  axes 
respectively.  Initial  values  of  -15  fig,  -150  fig,  and  +125  fig  respectively,  ap¬ 
proached  the  values  in  Figure  5  in  about  one  day. 

(b)  The  approximately  2-ftg  peak-to-peak  noise  on  the  radial  and  cross-track 
axes  precluded  obtaining  useful  data  from  them.  The  noise  in  these  axes  was  not 
satisfactorily  explained.  However,  no  evidence  of  a  feature  corresponding  to  the 
Z-axis  spike  was  found. 

Several  orbital  maneuvers  were  performed  to  study  further  the  instrument 
behavior.  These  tests  showed  that  all  three  axes  were  properly  responding  to 
atmospheric  density  and  to  vehicle  dynamics.  However,  because  of  the  contamin¬ 
ation  of  the  data  by  the  temperature  anomaly,  it  was  determined  that  geophysical 
information  could  not  be  extracted  with  sufficient  accuracy  to  extend  present  model¬ 
ling  efforts. 

The  contractor  (Bell  Aerospace  Textron)  determined  that  the  high  initial  bias 
values  were  due  to  charge  build-up  on  the  electrode  carrier  dielectrics.  They  also 
suggested  that  increasing  the  sensitivity  would  result  in  lower  steady-state  bias 
values.  The  initial  scaling  had  been  limited  to  the  maximum  sensitivity  the  con¬ 
tractor  recommended.  They  now  felt,  based  on  flight  experience,  that  the  sensi¬ 
tivity  could  be  improved. 


The  cause  of  the  temperature  anomaly  was  accidently  discovered  during  tests 
of  the  next  flight  unit  at  Bell  Aerospace  Textron.  It  was  found  that  application  of 
heat  to  the  proof-mass  housing  caused  a  jump  in  the  Z-axis  output.  Flexing  of  the 
thin  housing  walls  in  response  to  temperature  variations  created  a  flow  of  air  out 
of  a  vent  hole  in  the  electrode  carrier  and  against  the  cylinder  flange  that  detects 
Z-axis  motions.  The  geometry  of  the  proof -mass  housing  area  and  the  air  flow 
are  shown  in  Figure  6.  This  diagram  also  shows  that  die  X  and  Y  axes  were  ap¬ 
parently  not  affeci'.i  since  the  flow  is  perpendicular  to  them.  On  the  single -axis 
instruments  the  vent  hole  was  placed  symmetrically  with  respect  to  die  flange  on 
the  stainless  steel  spacer  ring.  In  the  development  of  die  three -axis  system  die 
hole  could  not  be  drilled  through  the  thin  (0.  056  cm)  spacer  ring,  now  made  of 
ceramic,  hence  it  was  placed  on  the  electrode  carrier. 


SEALED  COVER 


Figure  6.  Proof-Mass  Housing  Geometry 


Based  on  flight  data  and  the  testing  at  Bell  Aerospace  Textron,  the  following 
changes  were  incorporated  into  the  next  flight  unit: 

(a)  A  second  vent  hole  was  drilled  into  die  electrode  carrier.  It  was  placed 
symmetrically  opposite  the  first  hole  with  respect  to  the  proof -mass  flange. 


(b)  A  temperature  controller  was  added  to  eliminate  the  need  for  a  temperature- 
dependent  scale  factor  correction  to  the  flight  data.  This  controller,  designed  to 
operate  at  105°F  ±  0. 3°F  consisted  of  two  1.  5-W  heaters.  A  section  was  added  at 
each  end  of  the  housing  to  accommodate  the  heaters. 

(c)  The  sensitivity  was  increased  for  all  three  axes  as  shown  in  the  data  of 
Table  2.  This  was  done  to  reduce  the  bias  as  well  as  to  achieve  more  accurate 
measurements. 

(d)  Low  noise  preamplifiers  were  used  to  reduce  noise  in  the  cross  axes. 

This  modified  instrument  was  operated  during  1979.  Early  flight  data  indicated 

that  the  revisions  put  into  this  unit  resulted  in  extremely  accurate  acceleration 
measurements  on  all  three  axes.  Results  are  shown  for  one  orbit  in  Figure  7. 

This  format  is  identical  to  that  of  Figure  5.  The  Z-axls  data  show  a  drag  profile 
with  maximum  value  at  minimum  altitude  and  a  7 -tig  amplitude.  There  is  no  evi¬ 
dence  of  the  spike  or  other  anomalous  temperature  effects.  The  Y-axis  noise  is 
sufficiently  low  that  die  orbit-attitude  pulses,  more  frequent  near  perigee,  are  dis¬ 
tinguishable.  On  the  X-axis  a  1-tig  peak- to -peak  modulation  is  due  to  atmospheric 
wind  motions.  Stabilization  of  the  instrument  operating  temperature  at  110°F 
±  0.  3°F  occurred  within  one  hour  after  instrument  turn-on.  Table  3  summarizes 
the  improvements  in  sensitivity,  noise  reduction,  and  bias  achieved  with  this  flight. 
Applications  of  these  high  accuracy,  high  resolution  measurements  are  discussed 
in  Section  5. 

Normal  operation  continued  for  about  22  days.  Stabilization  of  the  tempera¬ 
ture  5°F  above  the  value  obtained  in  ground  tests  was  attributed  by  Bell  Aerospace 
Textron  to  the  effect  of  convection  of  heat  away  from  the  proof -mass  housing  dur¬ 
ing  ground  tests.  Since  they  had  not  considered  tills  before  launch,  data  were  re¬ 
duced  by  using  a  linear  extrapolation  of  the  calibration  data. 

Instrument  performance  began  to  degrade  after  about  22  days  of  operation. 

The  problem  initially  showed  up  as  a  drift  in  the  Y-axis  bias.  The  other  two  axes 
were  affected  but  to  a  lesser  extent.  This  drift  increased  with  time  on  all  three 
axes.  Figure  8  shows  a  plot  of  the  response  of  all  three  axes  after  about  one  week 
of  anomalous  operation.  The  sudden  shifts  toward  a  zero  bias  were  frequently  cor¬ 
related  with  vehicle  noise.  This  problem  persisted  throughout  the  flight.  Events 
such  as  orbit-adjusts  and  instrument  turn-off/turn-ons  lowered  the  bias  and  its 
drift  rate. 

Analyses  of  characteristics  of  the  anomalous  Instrument  operations  were  pro¬ 
vided  to  tiie  contractor  but  they  were  unable  to  account  for  these  phenomena.  Tests 
of  various  electrical  components  were  performed  but  the  bias  drift  could  not  be 
reproduced.  A  possible  problem  discovered  during  testing  was  that  the  two  sealed 
covers  on  the  ends  of  the  proof-mass  housing  could  be  assuming  a  potential.  This 
would  occur  if  they  became  electrically  Isolated  from  the  case.  Tests  of  one 


RAW  SETA  ACCELERATIONS 


DAY  79094 


10 

5 

0 

-5 


-10 


-IS 


10 

5 

0 

-5 


-10 


-15 


72000  74000  76000  78000 


TIME  (SEC) 


Figure  7.  Raw  Data  for  One  Orbit,  Day  79094 
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Table  3.  Instrument  Improvement  Summary 
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pre-flight  unit  were  performed.  An  open  circuit  condition  was  found  between  the 
case  and  one  cover.  Hie  cause  was  attributed  to  failure  of  the  epoxy  seal  induced 
by  internal  pressure  or  continuous  flexing  due  to  heater  cycling.  If  this  occurred 
in  orbit  the  accelerometer  output  could  exhibit  the  observed  erratic  behavior. 
Another  potential  problem  found  was  that  Interference  occurred  between  die  voltage 
regulator  switching  frequency  at  various  line  voltages  and  the  proof-mass  suspen¬ 
sion  frequencies  (60, 80,  and  100  MHz  for  X,  Y,  and  Z  respectively).  This  inter¬ 
ference  results  in  a  step  change  in  output  with  no  subsequent  time  variations. 

Hence  this  possible  problem  is  distinct  from  the  bias  drift  anomaly. 
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Figure  8.  Bias  Drift  Data  for  One  Day.  Day  78108 

Since  the  teat  flight  did  not  exhibit  a  bias  problem.  poBBible  effects  of  the  four 
changes  incorporated  after  this  flight  were  examined.  The  second  vent  hole  and 
the  low  noise  pre -amplifiers  were  determined  to  be  unrelated  to  the  problem.  The 
increase  in  C-range  sensitivity  was  suspected  as  a  possible  cause  of  an  instability 
in  the  instrument  response  to  acceleration  impulses,  and  the  duty-cycle  heaters 
were  suspected  as  a  possible  cause  of  intermittent  electrical  interference.  Pre¬ 
flight  tests  of  the  instrument  did  not  reveal  any  anomalies.  However,  because  of 
the  limitations  of  testing  in  earth's  one-g  environment  and  the  possibility  of  a  com¬ 
ponent  degradation  in  orbit  leading  to  marginal  instrument  performance,  the  latter 
two  changes  are  not  assumed  to  be  unrelated  to  the  bias  drift  problem. 

As  a  result  of  the  above  tests  and  analyses  of  flight  data,  it  was  decided  that 
the  following  changes  would  be  implemented  into  any  future  flight  units: 

(a)  Ground  the  cover  to  the  case. 

(b)  Replace  the  switching  (voltage)  regulator  with  a  nonswitching  type. 

(c)  Use  the  same  sensitivity  levels  as  those  on  the  test  flight  on  all  three  axes. 

(d)  Disconnect  the  heaters  internally. 


4.  DATA  ANALYSIS 


The  accelerometer  output  contains  contributions  related  to  vehicle  dynamics 
and  instrument  operating  characteristics  as  well  as  those  due  to  aerodynamic  drag 
and  winds.  The  vehicle  dynamic  accelerations  (ay)  are  mainly  those  due  to  space¬ 
craft  vibrations  induced  by  the  attitude  control  system  (discussed  in  Section  2). 
Centripetal  acceleration  and  the  gravity  gradient  effect,  each  of  the  order  of  0.  32  ng 
but  oppositely  directed,  provide  a  negligible  (<0.  01  fig)  resultant  dynamic  acceler¬ 
ation  in  the  along-track  axis.  Also,  the  radial -axis  data  contain  essentially  con¬ 
stant  contribution  due  to  the  gravity  gradient  acceleration  (-0. 08  fig).  Extraneous 
accelerometer  outputs  consist  of  a  bias  (Sg)  and  random  noise  (a^ ).  With  the 
centripetal  and  gravity  gradient  accelerations  Incorporated  as  part  of  the  bias 
term,  the  average  accelerometer  output  over  the  2. 045-sec  sample  time  (a^,)  is 
then  represented  by: 

aTi  =  aDi  +  avi  +  aBl  +  aNl  •  (1) 

where  i  =  X,  Y,  Z  since  the  various  acceleration  sources  can  be  different  along 
each  axis.  The  aerodynamic  component  is  derived  using  Eq.  (1)  in  the  form: 

aDl  =  aTi  "  favi  +  aBi  +  aNl^  * 

Numerical  filtering  techniques  (for  example,  Noonan  et  al5)  are  used  to  remove 
the  high  frequency  components  of  a^  due  to  satellite  dynamics  (av)  and  Instrument 
noise  (aN).  Figure  9  shows  the  result  of  applying  a  filter  passing  only  frequencies 
less  than  0. 01  Hz  to  the  data  of  Figure  7.  The  data  are  offset  from  zero  since  fil¬ 
tering  does  not  remove  the  nearly  steady-state  bias.  Since  Die  in-track  acceler¬ 
ometer  axis  senses  aerodynamic  drag  throughout  the  orbit,  bias  cannot  be  directly 
deduced  by  the  Z-axis  accelerometer  output  at  apogee.  A  bias  value  is  estimated 
for  each  orbit  essentially  by  subtracting  a  drag  acceleration  based  on  a  model  den- 
sity  value  from  Die  filtered  accelerometer  output  at  maximum  altitude.  This 
density  typically  corresponds  to  a  drag  value  of  0. 4  fig.  For  the  cross-track  and 
radial  axes,  aerodynamic  accelerations  are  negligible  at  this  altitude.  However 
bias  determination  is  complicated  by  low  frequency  noise  (<0.  01  Hz)  with  an 

5.  Noonan,  J.  P„  Ftoretti,  R.W. ,  and  Hass,  B.  (1975)  Digital  Filtering  Analysis 

Applied  to  Ihe  Atmosphere  Explorer -C  Satellite  MESA  Accelerometer  Data. 

AFCRL- 13-95-0393;  ADA  0l5?65. - 

6.  Jacchla,  L.G.  (1971)  Revised  Static  Models  of  the  Thermosphere  and  Exo- 

sphere  With  Empirical  Temperature  Profiles.  Spec.  Rep.  332.  Smithson. 
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amplitude  of  about %  04  pg  to  the n#9  ttad  9.=08  pg  tn  tbe  X  attd  Y  axes.  Hie  - 
aohrce  of  the  hoi#*  appears  tb  be  satellite  vibration  tviteteemeasOTedapeetrum 
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estimate.  First,  a  SOOnteebiat  average  centeredoft  maximum  altitude  was  used 
Instead  of  ]«il  themaxtmom  altttadeYalae.  Second,  the  average  bias  estimates 
for  ail  data  daring  22  days  of  normal  operation  -were  III  with  an  expression  of  the 
form:  ^ 


a^ft)  ■  Kj  +  ICj  exp  (Kj  t)  +  exp  (K^t) 


where  %  » time  in  days  from  day  79077. 09,  and  the  coefficients  for  each  axis  are: 


K1 

-147.3229  X 10-6 

2.  5344  X  10"6 

4.4670  X 10" 

*2 

9. 1053  X 10"6 

-6.4650  X 10"6 

-1.4235  X  10" 

*3 

-1.6944 

-1.5259 

-0. 7675 

K4 

147.  7283  X  10"6 

0.4453  X 10"6 

-1.8088  X  10"' 

K5 

-0. 1463  X  10'4 

0.0583  X  10‘6 

-0. 1027 

Results  are  shown  in  Figure  10  for  each  axis.  The  solid-line  fit  to  the  individual 
points  represents  the ‘bias. 

The  equations  used  for  the  determination  of  atmospheric  density  and  wind 
velocities  are: 


aDt 


CjPV2 


(3) 


where 

i  -  X.  Y.  Z 
aDl  =  drag  acceleration 
Aref  =  satellite  frontal  area 
m  =  satellite  mass 


* 


C 


l 


*  atmospheric  density 


=  atmospheric  mass  velocity  relative  to  the  satellite 


(4) 


with  Vq,  V^,  and  Vw  representing  respectively,  inertial  satellite 
velocity,  the  atmospheric  rotation  velocity  (assumed  equal  to  earth's 
rotation  velocity)  and  the  neutral  wind  velocity 


=  drag  coefficients:  =  Cj  (Vx,  Vy.  Vz). 


7 

The  results  of  Sentman  are  used  to  describe  the  components  of  the  satellite  drag 
coefficient.  In  this  formulation  the  equation  for  the  component  in  a  particular 
direction  of  the  total  force  on  an  element  of  area  is  given  by: 


7.  Sentman,  L.H.  (1961)  Free  Molecular  Flow  Theory  and  Its  Application  to  the 
Determination  of  Aerodynamic  Forces.  LMSC-448514,  Lockheed  Missiles 
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Figure  10.  X-.Y-,  and  Z-axis  Bias  Values.  Days  79079-79100  (Cont) 


|c  -  |(5k  +  7i +nt)  [t  •  (1+erf  7S)+^=  e"72  s2J 

+  (1  +  erftS)  +  (1  +  erf  YS)  +  -4re_'v2s2]  j  dA  .  (5) 


where 


A„,PV4 


aerodynamic  force 


I,  y,  rj  -  direction  cosines  between  mass  velocity  vector  and  axes  of 
element  of  area,  respectively 


k,  i.t  =  direction  cosines  between  the  axes  of  element  of  area,  and  the 
direction  in  which  force  is  desired  (accelerometer  axis) 


erf  yS  *  error  function  of  yS 


=  molecular  speed  ratio  =  (mass  velocity  of  gas)/(most  probable 
random  speed  of  molecules) 


Tj  =  temperature  of  incident  molecules 

Tr  =  temperature  of  reflected  molecules  (assumed  -  300°K) 

Eq.  (5)  is  integrated  over  the  entire  satellite  surface  area  to  obtain  the  drag 
coefficient  for  a  particular  axis.  Values  appropriate  to  this  satellite  were  used 

g 

with  Eq.  (5)  to  derive  drag  coefficients.  Figure  11  shows  an  example  of  drag 
coefficients  for  each  axis  as  a  function  of  pitch  angle.  These  results  were  ob¬ 
tained  by  assuming  a  temperature  (T^  of  1000°K  and  a  total  satellite  velocity  of 
8  km/sec  along  the  Z  axis.  The  values  have  been  normalized  by  the  satellite 
Z-axls  cross-sectional  area.  Negative  values  Indicate  that  the  drag  force  is 
along  the  direction  defined  as  negative  for  that  axis. 
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Figure  11.  Drag  Coefficients  for  X,  Y,  and  Z  Directions 


In  the  solution  of  Eq.  (3),  only  three  acceleration  measurements  are  available 
(a{)  while  there  are  four  unknowns  (p,  V).  An  iterative  technique  has  been  develop¬ 
ed  to  determine  density  and  winds.  The  Z-axls  satellite  velocity  (~8  km/sec)  is 

8.  Mazzella,  A.J.,  Jr.,  (1982)  Drag  Coefficients  for  a  Satellite.  A FOL  Technical 
Report  (in  press). 


large  compared  to  its  uncertainty,  hence  an  initial  estimate  of  the  density  is  ob¬ 
tained  by  assuming  VQ  ■  -VG  +  VA  and  calculating 


» 0 


DZ 


V?  C2(V0) 


(6) 


This  density  is  used  to  calculate  the  three  velocity  components,  V,  and  to  obtain 
the  wind  as  Vw  =  V  -  VQ.  To  refine  the  solution  a  new  V  =  fV^,  Vy,  Vq2)  is  used 
to  calculate  a  new  p.  This  new  p  is  then  used  to  calculate  new  velocities.  The 
process  is  repeated  until  the  value  of  18  le8S  than  0.  02  |VZI  • 

The  errors  in  the  calculated  values  of  density  are  attributed  to  accelerometer 
bias,  attitude  uncertainties,  drag  coefficient  estimates,  noise  within  the  desired 
signal  frequency  spectrum,  and  the  inseparability  of  density  and  winds  along  the  in¬ 
track  axis.  Systematic  errors  due  to  the  drag  coefficient  are  estimated  to  be  within 
±5  percent,  while  those  due  to  bias  vary  from  ±0.4  percent  at  170  km  to  ±  2  per¬ 
cent  at  220  km.  Low  frequency  noise  (0. 04  Mg)  contributes  a  random  error  of 
±  0.  5  percent  at  170  km  and  ±2.5  percent  at  220  km.  In  addition  an  assumed  wind 
uncertainty  of  200  m/sec  along  the  Z  axis  contributes  a  random  error  of  ±  5  per¬ 
cent.  Attitude  errors  are  negligible.  For  applications  of  the  data  requiring  only 
knowledge  of  the  drag,  rather  than  the  density,  the  errors  are  further  reduced. 

The  errors  due  to  drag  coefficient  and  Z-axis  winds  are  eliminated.  Further,  if 
long-term  average  drag  (order  of  one  orbit)  is  required  the  low  frequency  noise 
average  will  be  approximately  zero.  The  major  error  contribution  to  cross-axis 
wind  determination  is  the  low  frequency  noise.  The  amplitude  of  this  noise  is 
0.  08  ng  while  the  maximum  acceleration  due  to  wind  is  less  than  0. 1  Mg.  As  a 
result,  there  is  a  random  error  of  ±  60  m/sec  at  170  km  and  ±  160  m/sec  at 
200  km.  Uncertainties  in  attitude  and  the  Z-axis  wind  result  in  an  additional 
random  error  of  about  ±  10  m/sec  at  all  altitudes.  By  combining  results  from 
n  orbits  obtained  during  similar  solar  and  geophysical  conditions  an  average  wind 

behavior  can  be  derived  with  the  error  reduced  by  a  factor  of  n"1^2.  A  complete 

q 

summary  of  the  density  and  wind  error  analysis  has  been  given  elsewhere. 


5.  GEOPHYSICAL  STUDIES 

Global  studies  of  present  and  future  SETA  density  and  wind  data  will  lead  to 
an  improved  understanding  of  atmospheric  dynamics  and  consequently  to  more 
accurate  models.  Upper  atmosphere  heating  above  150  km  is  mainly  by: 

9.  Mazzella,  A.J. ,  Jr.  (1982)  Error  Analysis  for  Densities  and  Winds,  AFGL 
Technical  Report  (in  pressT 


(1)  solar  electromagnetic  energy  at  extreme  ultraviolet  wavelengths, 

(2)  energy  from  the  solar  wind  plasma  resulting  in  particle  heating  and  convec¬ 
tion  electric  fields  (Joule  heating),  and 

(3)  upward  propagating  gravity  waves  and  tides. 

Gravity  waves  and  winds  are  involved  in  transporting  energy  from  the  auroral  zone 
to  lower  latitudes.  The  SETA  data  permit  investigation  of  aspects  of  each  of  these 
effects.  Section  5. 1  gives  a  brief  overview  of  these  heat  sources  and  describes 
how  the  satellite  orbital  characteristics  relate  to  their  study.  Sample  measure¬ 
ments  from  the  operational  SETA  flight  are  applied  to  problems  involving  each 
heat  source  in  turn  in  Section  5.2.  Section  5. 3  demonstrates  the  capability  of  the 
accelerometer  to  measure  neutral  winds.  Although  treated  separately  here,  simul¬ 
taneous  wind  and  density  data  obtained  on  an  every  orbit  basis  over  a  wide  latitude 
region  for  the  first  time,  will  be  utilized  in  planned  geophysical  studies. 

S.1  Relation  of  Data  to  Atmospheric  Studio 

5.  1.  1  NEUTRAL  ATMOSPHERIC  HEATING  AND 
DYNAMICS  OVERVIEW 

The  main  source  of  thermospheric  heating  Is  the  solar  electromagnetic  energy. 
It  has  major  variations  related  to  the  eleven-year  solar  cycle  and  the  twenty-seven 
day  solar  rotation  period.  This  energy  is  deposited  mainly  at  low  latitudes  and  in 
the  summer  hemisphere.  The  density  maximizes  in  the  vicinity  of  the  sub-solar 
point  and  reaches  a  minimum  value  near  the  anti-solar  point.  This  diurnal  vari¬ 
ation  is  modulated  by  Die  seasonal  migration  of  the  sun's  declination.  Jacchia6 
contains  an  excellent  description  of  these  variations.  This  density  structure  pro¬ 
duces  pressure  gradients  that  cause  an  atmospheric  circulation.  An  early  wind 
model  based  on  pressure  gradients  was  developed  by  Kohl  and  King.  10  A  wind  flow 
is  set  up  essentially  from  the  dayside  to  the  nightatde  of  the  earth.  Near  the  equa¬ 
tor  the  winds  are  mainly  zonal  (eastward  for  local  times  later  than  the  time  of  the 
diurnal  maximum,  westward  for  earlier  local  times).  At  higher  latitudes  the  wind 
is  predominately  toward  the  pole  on  the  dayside  and  away  from  the  pole  on  the 
nights ide.  Typical  wind  speeds  associated  with  the  solar  EUV  heat  source  are  of 
the  order  of  100  m/sec. 

The  solar  wind,  in  contrast  to  the  solar  EUV,  provides  a  heat  source  that  is 
effective  at  high  latitudes.  The  earth's  approximately  dipole  magnetic  field  is 
radically  modified  beyond  a  few  earth  radii  through  interaction  with  the  solar  wind 
plasma  that  flows  by  earth  with  a  velocity  in  excess  of  400  km/sec.  The  effect  of 
this  supersonic  flow  on  earth's  geomagnetic  field  is  to  compress  it  on  the  dayside 

10.  Kohl,  H. ,  and  King,  J.  W.  (1967)  Atmospheric  winds  between  100  and  700  km 
and  their  effects  on  the  ionosphere.  J.  Atmos.  Terr.  Phys.  29:1045. 


and  elongate  it  (beyond  the  moon's  orbit)  on  the  nightside.  This  deflection  of  the 
magnetic  field  creates  a  cavity  to  which  the  solar  wind  does  not  have  direct  access. 
This  cavity  and  its  boundary  regions  are  known  as  earth's  magnetosphere.  A 
schematic  view  of  the  magnetosphere  is  shown  in  Figure  12. 11  Its  interior  con¬ 
sists  of  several  regions  containing  plasma  with  a  variety  of  particle  densities  and 
energies.  Electrodynamic  coupling,  or  interaction,  occurs  between  the  magneto- 
sphere,  ionosphere,  and  upper  atm  osphere.  The  mechanisms  are  not  well  under¬ 
stood.  It  is  known  that  the  magnetosphere  affects  the  neutral  upper  atmosphere 

mainly  by  driving  large-scale  currents  through  the  ionosphere.  These  currents 
13 

result  in  Joule  heating  at  auroral  zone  latitudes.  Both  the  magnitude  and  spatial 
extent  of  the  high  latitude  heating  correlate  with  geomagnetic  activity  and  are  ex¬ 
tremely  variable.  As  a  result,  large  and  sudden  density  variations  occur.  Solar 
wind  related  energy  is  also  deposited  directly  by  particle  precipitation.  One 
mechanism  involves  electrons  that  are  accelerated  in  the  geomagnetic  tail  and 
precipitate  at  high  latitudes.  The  polar  "cusp"  region  (top  of  Figure  12)  provides 
another  means  of  particle  precipitation.  The  cusp  appears  to  be  like  a  hole  in  the 
magnetosphere,  permitting  solar  wind  particles  to  directly  enter  the  upper  atmos¬ 
phere.  The  affected  region  corresponds  roughly  to  the  noontime  auroral  oval. 

Strong  wind  systems  can  be  induced  by  the  high  latitude  heating  mechanisms. 
The  circulation  is  driven  from  high  latitudes  to  low  latitudes.  Its  strength 

depends  on  the  magnitude  of  the  heat  Source.  During  large  storms  this  circulation 

14  15 

can  reverse  that  due  to  solar  EUV  heating  (for  example,  Straus  et  al,  Straus, 

1 6 

Dickinson  et  al  ).  These  winds  are  a  major  mechanism  for  heat  transport. 

Atmospheric  tides  are  regular  large-scale  oscillations  excited  by  gravitational 
forces  of  the  sun  and  moon  or  by  absorption  of  solar  radiation.  They  have  a  basic 
period  of  one  lunar  or  one  solar  day  with  higher  harmonics  down  to  a  period  of  six 
lunar  or  solar  hours.  Previous  accelerometer  data  have  shown  that  below  200  km 
the  amplitude  of  the  solar  semidiurnal  variation  exceeds  that  of  the  diurnal 
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Figure  12.  Diagram  of  Earth's  Magnetosphere 
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bulge  *  and  that  toe  phase  of  toe  diurnal  variation  shifts  to  later  local  times  from 

18 

the  equator  to  mid-latitudes.  Tidal  phenomena  have  been  reviewed  recently  by 
Forbes  and  Garrett.18 

Large-scale  gravity  waves,  having  periods  different  from  those  of  tidal  waves 
and  not  directly  related  to  the  earth -moon-sun  geometry,  also  occur  in  toe  atmos¬ 
phere.  These  waves  may  be  generated  either  locally  in  the  thermosphere  or  near 
the  earth's  surface.  Those  originating  in  toe  thermosphere  are  frequently  associ¬ 
ated  with  toe  energy  released  in  a  narrow  zone  at  auroral  latitudes  during  geomag¬ 
netic  storms.  Theoretical  studies  have  shown  that  the  atmosphere  acts  as  a  filter, 

allowing  only  longer  wavelengths  (hundreds  of  km)  to  propagate  a  significant  dis- 

20  2 1 

tance  from  toe  auroral  zone.  *  Waves  propagating  into  toe  thermosphere  from 

below  have  been  postulated  to  result  from  earth's  orographic  features,  weather 

fronts,  and  let  stream  turbulence.  Evidence  of  a  tropospheric  source  for  wavelike 

22 

features  in  density  profiles  has  been  presented  by  Rice  and  Sharp. 

The  numerous  heating  and  circulation  processes  involved  in  determining  at- 

23 

mospheric  variability  are  summarized  in  Figure  13.  (This  figure  also  includes 

Because  of  toe  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  pages  51  and  52. 
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Figure  13.  Heat  Sources  tor  the  Upper  Atmosphere 


processes  outside  the  scope  of  the  present  discussion).  Both  neutral  and  ionized 
constituents  exhibit  regular  global  patterns  as  well  as  significant  short-term 
variations.  During  geomagnetically  quiet  times  the  thermosphere  is  driven  mainly 
by  solar  electromagnetic  energy.  High  latitude  processes  dominate  during  large 
geomagnetic  storms.  At  other  levels  of  geomagnetic  activity  an  intermediate  con¬ 
dition  prevails.  To  develop  more  accurate  models  our  knowledge  of  the  magnitude 
and  distribution  of  atmospheric  heating  must  be  improved  considerably. 
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5.  1. 2  ORBITAL  DATA  COLLECTION 

The  normal  data  collection  mode  for  die  accelerometer  is  continuous  operation 
in  die  most  sensitive  (C)  range.  This  permits  drag  acceleration  measurements  to 
be  obtained  over  a  complete  orbit  on  every  orbit.  High  accuracy  density  data  are 
obtained  over  about  a  15  o'* -latitude  band,  and  winds  are  obtained  over  about  a 
100°-latitude  band  with  an  accuracy  of  the  data  that  is  a  function  of  satellite  altitude, 
as  discussed  in  Section  4.  Figure  14  summarizes  die  orbital  characteristics  with 
respect  to  die  northern  hemisphere  auroral  zone  for  a  one  day  period.  Geomag¬ 
netic  latitude  and  magnetic  local  time  coordinates  are  used,  and  are  indicated  on 


Figure  14.  Intersection  of  Satellite  Tracks  for  One  Day  With  Northern  Hemi¬ 
sphere  Auroral  Zone 


34 


the  figure.  This  permits  plotting  the  satellite  track  with  respect  to  a  fixed  auroral 
zone.  Because  of  the  high  satellite  inclination  and  relatively  low  eccentricity,  both 
the  day  and  night  auroral  regions  are  sampled  on  every  orbit  at  approximately  the 
same  altitude. 

The  combination  of  continuous  operation  for  an  extended  period  of  time  over  a 
wide  latitude  region  with  routine  high  temporal  resolution  auroral  zone  data  has  not 
been  achieved  with  previous  accelerometer  flights.  These  factors  provide  a  unique 
opportunity  to  investigate  the  role  of  the  various  heat  sources  on  the  neutral 
atmosphere. 

5.2  Atnoapheric  Denrity  RmuUi 

An  example  of  an  atmospheric  density  profile  is  shown  in  Figure  15.  Density 
“3 

(g  cm  )  is  plotted  as  a  function  of  the  following  parameters:  GMT  (seconds), 
geographic  latitude,  geomagnetic  latitude,  local  time  (hours  and  minutes).  GMT 
(hours  and  minutes),  and  altitude  (km).  These  results  were  obtained  during  the 
geomagnetic  storm  to  be  discussed  in  Section  5. 2.  2.  Superimposed  on  the  general 
increase  of  density  with  decreasing  altitude  are  significant  oscillations.  These 
wave  structures  are  described  farther  in  Section  5.2.3. 
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Figure  15.  Atmospheric  Density  as  a  Function  of  Altitude,  Latitude,  and  Local  Time 


5.2. 1  SOLAR  FLUX  EFFECTS 

Thermospheric  heating  by  solar  EUV  radiation  is  represented  by  a  strong 
positive  correlation  of  density  and  the  solar  radio  flux  measured  at  10.7  cm 
(F1Q>7)  in  all  current  atmospheric  models.  The  radio  flux  is  assumed  to  track 
variations  in  EUV  flux  and  thus  reflect  variations  in  the  thermospheric  energy 
input.  The  long-term  component  (primarily  related  to  the  11 -year  solar  cycle) 
is  represented  by  a  flux  averaged  over  several  (usually  three)  solar  rotation  peri¬ 
ods  (F).  The  short-term  component  related  to  dally  changes  In  solar  activity  is 
represented  by  (F-F)  where  F  has  a  lag  of  about  one  day.  An  empirical  relation¬ 
ship  between  density  and  solar  flux  has  been  derived  for  low  solar  flux  conditions 
(F10. 7  <  100)  by  Marcos  and  Champion2*  using  regression  techniques  described 
in  Marcos  et  al.25  An  extension  to  high  solar  flux  conditions  based  on  measured 
data  is  required  since  a  linear  extrapolation  may  not  be  an  accurate  representation. 

Figure  16  shows  the  trend  of  the  density  variability  at  200  km  with  solar  flux 
as  determined  from  accelerometer  data.  The  solid  line  represents  die  relation  de- 
rived  for  F1Q  7  £  100  from  Marcos  and  Champion.  The  dashed  line  is  a  linear 
extrapolation  of  this  result.  Circles  represent  data  obtained  during  several  days 
in  April  1978,  a  period  of  moderate  solar  activity,  with  the  ROCA  experiment  on 

nil 

the  S3 -4  satellite.  Data  were  obtained  over  a  five -month  period,  but  due  to  pro¬ 
cessing  problems  at  SAMTEC,  delivery  of  data  tapes  to  AFGL  were  considerably 
delayed.  Routine  processing  of  the  data  commenced  during  CY  1981.  SETA  data 

for  several  days  in  April  1979.  a  period  of  high  solar  flux,  are  shown  as  crosses. 

27 

For  comparison,  values  predicted  by  the  MS  IS  model  are  shown  as  triangles. 

The  accelerometer  data  Indicate  a  larger  response  to  solar  flux  than  the  MSIS 
model.  A  quantitative  empirical  estimate  of  the  density  increase  per  unit  solar 
flux  increase  will  be  derived  with  additional  ROCA  and  SETA  data.  The  result 
can  be  directly  utilized  to  extend  the  present  empirical  model. 
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SOLAR  FLUX  ( xlOc*  W*tiT  •  Ha” 1 ) 

Figure  16.  Atmospheric  Density  as  a  Function  of  Solar 
Flux 


5.2.2  GEOMAGNETIC  ACTIVITY  EFFECTS 

Rapid  global  fluctuation  of  atmospheric  density  with  a  typical  duration  of  one 
to  several  days  are  connected  with  the  solar  wind  related  geomagnetic  disturbances. 
The  amplitude  of  the  density  increase  generally  depends  on  the  strength  of  the  geo¬ 
magnetic  storm  as  measured  by  the  three-hourly  planetary  index  Kp.  The  atmos¬ 
pheric  response  occurs  first  at  auroral  zone  latitudes  and  then  at  low  latitudes  with 
a  delay  of  about  6  to  8  h.  This  implies  that  the  energy  transferred  to  the  thermo¬ 
sphere  through  auroral  processes  at  high  latitudes  is  later  distributed  globally 
through  dynamic  processes.28  Recent  composition  measurements  and  theoretical 
studies  have  shown  that  atmospheric  winds  driven  by  the  high  latitude  heat  source 
produced  by  the  magnetic  storm  are  involved  in  the  global  redistribution  of  energy 
and  mass.28  The  phase  and  amplitude  of  the  atmospheric  response  as  a  function 
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of  latitude  are  not  well  modelled  at  present.  The  effect  also  varies  with  longitude, 
season,  local  time,  and  solar  flux. 

Figure  17a  shows  the  response  of  the  atmosphere  to  a  large  geomagnetic  dis¬ 
turbance.  Density  measured  at  80°,  60°,  30°,  and  0°  north  geographic  latitude 
at  about  1030  h  LT  and  Kp  are  plotted  as  a  function  of  time  for  the  period  1200  h, 
day  79093  to  2400  h,  day  79094  GMT.  The  curves  therefore  show  density  values 
obtained  on  the  dayside  of  the  atmosphere  during  24  consecutive  orbits.  The  alti¬ 
tude  at  a  given  latitude  varies  slightly  during  this  period  due  to  perigee  precession. 
The  densities  at  each  latitude  have  been  normalized  to  a  constant  altitude,  indicated 
for  each  curve. 

An  increase  in  Kp  starts  at  1500  h  on  day  93  and  rises  to  a  maximum  value  of 
8  for  the  three -hour  interval  beginning  at  2 100  h.  A  density  increase  occurs  at 
80°  and  60°  in  response  to  the  first  Kp  increase,  but  not  at  the  lower  latitudes. 

In  addition  to  this  phase  difference,  another  feature  is  the  greater  intensity  of  the 
thermospheric  response  at  high  latitudes.  At  lower  latitudes  the  density  continues 
to  increase  while  Kp  decreases.  A  mechanism  for  maintaining  these  enhanced 
densities  is  the  large-scale  circulation  pattern  with  rising  motions  at  high  latitudes 
and  descending  motions  at  middle  and  low  latitudes  described  by  Mayr  et  al.®° 

3 1 

Figure  17b  shows  the  atmospheric  response  predicted  by  the  Jacchia  1977 
model  that  utilizes  a  latitude -dependent  time  delay  for  the  geomagnetic  activity 
effect.  Model  values  at  the  same  latitudes  and  altitudes  as  Figure  17a  are  given. 

At  each  latitude  the  amplitude  and  phase  of  the  model  densities  follow  the  Kp  vari¬ 
ation  with  a  time  delay.  This  is  in  contrast  to  the  complex  density  structure  actu¬ 
ally  measured  at  high  latitudes  and  the  gradual  density  increase  found  at  lower 
latitudes.  Hence  the  model  gives  a  poor  representation  of  the  atmospheric  density 
response  to  this  geomagnetic  storm.  A  comprehensive  analysis  of  geomagnetic 
activity  variations  detected  with  SETA  data  has  been  initiated.  An  improved  de¬ 
scription  of  this  phenomenon  is  an  important  requirement  for  improved  atmospheric 
models. 

5.2.3  WAVE  STRUCTURES 

Gravity  waves  with  largest  amplitudes  observed  at  high  latitudes  were  detected 
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by  the  accelerometers  on  the  OV 1-15  satellite.  These  were  related  to  die 
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Figure  17a.  Atmospheric  Density  Response  to  a  Geomagnetic  Storm 
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Figure  17b.  Model  Prediction  of  Density  Response  to  a  Geomagnetic  Storm 


intensity  of  heating  as  indicated  by  the  auroral  electrojet  (AE)  index  by  Forbes  and 
33 

Marcos.  Theoretical  models  describing  the  spectrum  of  waves  generated  by 

20 

high  latitude  heating  show  occurrence  of  long  wavelengths  of  the  order  of  3000  km. 

Evidence  of  wave  structures  with  location  dependent  on  earth's  tropospheric  weather 

22 

patterns  have  been  found  in  Atmosphere  Explorer  satellite  data  by  Rice  and  Sharp. 
SETA  results  provide  a  unique  opportunity  to  statistically  distinguish  between  waves 
originating  in  the  lower  and  upper  atmosphere.  Their  respective  characteristics 
can  be  studied  over  a  wide  latitude  region  on  an  every-orbit  basis. 

Figure  18  shows  the  nature  of  wavelike  disturbances  associated  with  geomag¬ 
netic  activity.  Density  data  obtained  during  one  orbit  of  the  storm  period  discussed 
in  the  previous  section  are  plotted  as  a  function  of  time  and  latitude.  To  indicate 
better  the  wave  structure  features,  the  exponential  dependence  of  density  with  al¬ 
titude  has  been  removed  by  normalizing  the  data  to  an  altitude  of  200  km.  This 
facilitates  detection  of  the  wave  half -amplitude,  determined  as  the  deviation  from 
a  smooth  curve  obtained  through  the  density  data.  The  largest  half -amplitudes  of 
about  -22  percent  and  +27  percent  occur  at  15625  sec  GMT  and  16050  sec  GMT 
respectively.  In  geomagnetic  coordinates  these  maxima  both  correspond  to  about 
72°N  on  the  nightside  and  dayside  respectively.  The  data  clearly  show  a  decrease 
in  the  amplitude  of  the  oscillatory  behavior  with  latitude  away  from  the  auroral 
zone.  The  apparent  horizontal  wavelengths  are  estimated  by  the  time  between  suc¬ 
cessive  density  maxima  or  minima  multiplied  by  the  satellite  velocity.  Values  of 
the  order  of  1500  to  4000  km  are  observed.  It  is  concluded  that  these  results  are 
consistent  with  gravity  waves  generated  by  auroral  sources  during  a  period  of  high 
geomagnetic  activity.  This  sample  of  data  demonstrates  the  capability  of  the  accel¬ 
erometer  to  provide  knowledge  of  the  detailed  structural  atmospheric  features. 

These  measurements  are  required  for  a  more  comprehensive  understanding  of  at¬ 
mospheric  perturbations. 

5.3  Neutral  Atmospheric  Wind  Reaulta 

The  importance  of  neutral  winds  in  determining  the  global  density  structure 
was  Indicated  in  Section  5. 1.  There  have  been  few  satellite  measurements  of 
neutral  winds.  A  review  of  several  techniques  and  their  limitations  has  been  given 
by  Hanson  and  Heelis.  Samples  of  wind  data  were  obtained  from  experiments  on 
the  AE-C  satellite.  Results  reported  to  date  have  been  extremely  limited.  The 


33.  Forbes,  J.M. ,  and  Marcos,  F.A.  (1873)  Thermospheric  density  variations 

associated  with  auroral  electrojet  activity,  J.  Geophys.  Res.  78:3841. 

34.  Hanson,  W.B. ,  and  Heelis,  R.A.  (1975)  Techniques  for  me  .  ring  bulk  gas 

motions  from  satellites,  Space  Set.  Instrum.  1:493. 


GMT  19000  19290  19900  15790  16000  16290  16900  16790  17000  17290  17500 


6GLAT 

23 

40 

96 

73 

83 

71 

99 

37 

20 

3 

-14 

GMLAT 

34 

90 

69 

78 

74 

59 

44 

28 

9 

-7 

-23 

LT 

2203 

2132 

2133 

2046 

1922 

1128 

1050 

1033 

1023 

1015 

1007 

Figure  18.  Wave  Structure  in  an  Atmospheric  Density  Profile 
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Neutral  Atmosphere  Temperature  Experiment  has  detected  vertical  winds.  A 
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technique  for  obtaining  along-track  winds  with  the  Open  Source  Spectrometer  was 
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described  by  Knutson  et  al.  Cross -axis  winds  from  the  MESA  accelerometer 
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have  been  derived  for  several  orbits  by  Hanson  et  al.  Accelerometer  horizontal 
wind  measurements  were  also  obtained  during  part  of  the  four  day  flight  of  the 
LOGACS  experiment.40' 4 1  The  SETA  accelerometer  has  demonstrated  a  capability 
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to  measure  cross  -axis  winds  over  a  wide  latitude  region.  Analysis  of  these  results 
will  significantly  enhance  our  understanding  of  global  transport  processes  and  con¬ 
sequently  atmospheric  density  variability.  The  wind  data  are  discussed  separately 
in  this  section  to  emphasize  their  significance.  However,  these  results  will  be 
most  fully  exploited  by  analyzing  them  in  conjunction  with  the  simultaneously 
measured  density  data. 

Two  wind  profiles  derived  from  cross-track  data  are  given  in  this  section. 

These  horizontal  winds  were  measured  near  the  end  and  at  the  maximum,  respec¬ 
tively,  of  the  geomagnetic  storm  on  day  7  9094.  The  major  features  of  the  observed 
winds  are  described.  More  comprehensive  analyses  will  be  provided  in  forthcom¬ 
ing  reports.  Vertical  acceleration  data  obtained  for  the  same  time  periods  are 
also  shown.  At  present,  these  have  not  been  reduced  to  wind  values.  A  detailed 
study  of  possible  satellite  motions  in  the  radial  direction  is  being  made  to  ensure 
extraction  of  vertical  atmospheric  motions  with  high  accuracy.  Vertical  wind  re¬ 
sults  will  also  be  the  subject  of  a  later  report. 

Because  the  orbit  is  not  strictly  polar,  the  angle  between  the  zcnal  wind  com¬ 
ponent  and  the  accelerometer  sensing  axis  is  a  function  of  geographic  latitude. 
Figure  19  shows  this  relationship.  As  die  pole  is  approached,  the  component  of 
cross-axis  wind  measured  shifts  from  mainly  zonal  toward  meridional.  Figure  20 
shows  filtered  acceleration  data  from  the  X  axis.  A  three -orbit  period  on  day 
79094,  when  Kp  was  moderately  low,  is  presented  to  demonstrate  the  repeatability 
of  these  data.  The  rotation  velocity,  assumed  to  follow  that  of  die  earth,  is  about 
460  m/sec  at  the  equator  and  decreases  with  latitude  to  zero  m/sec  near  the  poles. 
(To  derive  the  horizontal  wind  component  the  accelerations  due  to  atmospheric 
rotation  are  removed  as  discussed  in  Section  4. )  The  aerodynamic  side  acceler¬ 
ations  induced  by  this  motion  can  be  seen  in  each  of  the  three  orbits.  The  maximum 
(negative)  acceleration  caused  by  the  eastward  atmospheric  motion  occurs  at  mid¬ 
latitude  near  minimum  altitude.  At  hitler  latitudes  the  amplitude  of  this  modula¬ 
tion  decreases.  This  is  due  to  the  slower  rotation  speed,  the  changing  angle  of 
attack  and  a  lower  density  at  the  corresponding  higher  altitudes.  As  the  satellite 
moves  over  the  pole  the  winds  due  to  atmospheric  rotation  are  directed  at  its 
opposite  side.  As  a  result  a  positive  acceleration  output  is  sensed. 

Neutral  horizontal  winds  derived  from  the  first  orbit  of  the  data  of  Figure  20 
are  shown  in  Figure  2 1.  Velocity,  in  units  of  m/sec,  is  plotted  as  a  function  of 
GMT,  geographic  and  geomagnetic  latitude,  local  time,  and  altitude.  Positive 
velocities  indicate  winds  blowing  from  west-to-east  on  the  dayslde  of  the  earth. 

On  the  nightslde  positive  velocities  indicate  westward  winds.  At  low  and  middle 
latitudes  the  daytime  winds  are  observed  to  be  mainly  west-to-east  with  velocities 
of  about  100  to  200  m/sec.  If  the  diurnal  bulge  were  located  at  1600  h,  the  pressure 
gradients  at  low  latitudes  and  a  longitude  corresponding  to  about  1030  h  would  drive 
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Figure  21.  Horizontal  Wind  Velocities  During  a  Period  of  Moderately  Low  Geomag¬ 
netic  Activity 

a  wind  in  the  opposite  direction  (east-to-west).  Hence  these  observed  winds  are 
not  consistent  with  a  sub-solar  density  bulge.  However,  as  noted  in  Section  5. 1. 
the  semidiurnal  tide  predominates  below  200  km  and  is  a  possible  mechanism  for 
the  observed  wind  direction.  At  higher  latitudes  from  about  70  to  83°gg  (geo¬ 
graphic  latitude)  on  the  dayside.  the  winds  are  generally  westward  with  speeds  of 
0  to  50  m/sec.  On  the  nightside  the  winds  are  again  generally  west-to-east.  How¬ 
ever,  these  results  are  dominated  by  a  shift  to  an  east-to-west  flow  between 
64. 8°gg  (68. 8  gm)  and  7 1.  5°gg  (75. 6°  gm)  with  a  maximum  velocity  of  277  m/sec 
at  68. 8°gg  (72. 9°gm).  These  large  wind  values  may  be  an  effect  of  the  diminishing 
geomagnetic  storm  rather  than  a  persistent  feature.  This  will  be  resolved  with  the 
reduction  of  more  data. 

Figure  22  shows  cross-axis  acceleration  data  obtained  earlier  on  day  79094 
for  a  three -orbit  period  when  the  level  of  geomagnetic  activity  was  high.  Consid¬ 
erable  structure  is  seen  at  high  latitudes  in  contrast  to  the  data  of  Figure  20.  The 
wind  profile  derived  for  one  orbit  is  shown  in  Figure  23.  This  orbit  corresponds 
to  a  time  when  Kp  »  8,  the  maximum  value  for  this  storm.  Significant  structure  is 
observed  at  high  latitudes.  The  maximum  wind  velocity  for  this  orbit  is  597  m/ sec. 
It  occurs  on  the  nightside  of  the  earth  at  81. 6°gg  (78.  0°gm).  On  the  dayside  there 
is  a  sharp  transition  with  die  winds  shifting  from  westward  to  eastward  at  <77. 2°gg 
(65. 7°gm).  This  feature  was  also  observed  in  the  LOGACS  data  during  extremely 
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high  geomagnetic  conditions.  '  The  maximum  eastward  wind  is  475  m/sec  at 
72. 2°gg  (60. 7°gm).  At  latitudes  below  60°N  gg  the  wind  is  eastward  at  about 
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100  to  200  m/sec,  implying  that  the  storm  effects  are  confined  to  high  latitudes. 
Hernandez  and  Roble  have  suggested  that  a  large-scale  electric  field  associated 
with  magnetospheric  convection  could  explain  an  abrupt  transition  in  zonal  winds. 
The  observed  features  appear  to  be  related  to  large  amounts  of  energy  being  de¬ 
posited  locally  at  high  latitudes  and  being  subsequently  distributed  globally.  Pres¬ 
sure  gradients  associated  with  high  latitude  heating  significantly  modify  the  flow 
expected  from  solar  EUV  heating.  The  two  examples  of  horizontal  winds  (Fig¬ 
ures  2 1  and  23)  demonstrate  the  complexity  of  the  global  wind  field.  Analysis  of 
the  full  set  of  SETA  measurements  will  permit  resolution  of  many  features  of  the 
flow  of  energy  and  momentum  in  the  thermosphere. 

Figure  24  shows  radial  axis  filtered  acceleration  data.  Figure  24a  corre¬ 
sponds  to  the  time  period  of  Figure  20  and  Figure  24b  corresponds  to  the  time 
period  of  Figure  22.  Data  for  all  six  orbits  show  characteristic  "bumps"  occurring 
near  perigee  and  near  50°S.  These  bumps  correspond  to  an  increased  downward 
acceleration.  If  interpreted  as  vertical  velocities,  the  accelerations  of  Figure  24 
are  of  the  order  of  50  to  100  m/sec.  Presently,  vertical  speeds  of  the  order  of  a 
few  m/sec  have  been  predicted.  There  does  not  appear  to  be  a  significant  depend¬ 
ence  on  geomagnetic  activity.  No  attempt  to  interpret  these  results  will  be  given 
here.  The  data  are  of  extreme  interest  and  will  be  thoroughly  analyzed. 
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Geomagnetic  Activity 


6.  CONCLUSIONS 

An  accelerometer  has  been'  developed  using  a  single  proof  mass  to  measure 
simultaneously  satellite  accelerations  in  three  orthogonal  axes.  The  instrument 
is  based  on  the  design  of  the  highly  successful  single-axis  MESA  accelerometer. 
Future  instruments  will  be  modified  to  eliminate  bias  drift  anomalies. 

Specific  applications  of  the  extremely  accurate  data  obtained  during  a  22-day 
period  of  one  flight  demonstrated  the  importance  of  this  data  in  solving  geophysical 
problems.  Density  variations  related  to  long-term  (solar  cycle)  and  short-term 
(geomagnetic  activity  and  related  wave  propagation)  phenomena  were  shown  in 
Section  5  to  provide  valuable  inputs  to  improve  current  empirical  atmospheric 
models.  A  most  significant  development  of  this  flight  was  the  capability  to  detect 
wind  motion  for  the  first  time  over  a  wide  latitude  region  on  an  every-orbit  basis. 
Detailed  knowledge  of  die  spatial  and  temporal  wind  features  permit  extension  of 
theories  of  atmospheric  density  dynamics  and  will  provide  fundamental  knowledge 
of  the  magnetosphere -ionosphere -neutral  atmosphere  coupling.  Maximum  advan¬ 
tage  of  the  data  will  be  taken  by  performing  studies  of  the  simultaneous  features  of 
density  variability,  die  occurrence  of  wave  structure  and  the  neutral  wind  motions. 
Where  possible,  other  data,  such  as  Incoherent  scatter  (for  example,  Evans45) 
measurements  of  temperature  and  plasma  motions  will  be  incorporated  into  coordi¬ 
nated  analyses. 

45.  Evans,  J.  V.  (1975)  A  review  of  F  region  dynamics.  Rev.  Geophvs.  Space 
Phys.  13:887. 


Studies  of  data  from  the  SETA  flight  will  be  extended  to  examine  numerous 
other  atmospheric  structure  features.  These  include  horizontal  wind  and  neutral 
density  variations  as  a  function  of  latitude,  longitude,  and  local  time.  Data  ob¬ 
tained  in  die  high  latitude  regions  will  be  of  particular  importance  for  elucidation 
of  heating  characteristics.  The  accuracy  with  which  vertical  winds  can  be  extracted 
will  also  be  determined.  Their  measurements  would  provide  significant  new  infor¬ 
mation  regarding  the  circulation  of  die  upper  atmosphere. 

Atmospheric  dynamic  processes  have  been  shown  to  be  extremely  complex  and 
variable.  The  combination  of  accelerometer  capabilities  and  orbital  characteris¬ 
tics  provides  a  unique  opportunity  to  obtain  new  fundamental  data.  Long-term 
synoptic  measurements  anticipated  from  future  flights  coupled  with  the  present 
results  will  provide  a  major  contribution  toward  the  goal  of  developing  models  that 
can  accurately  predict  atmospheric  behavior. 
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